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Abstract 

We calculate the flux of cosmic positrons from the dark matter annihilation 
in the littlest Higgs model with T-parity. The dark matter annihilates mainly 
into weak gauge bosons in the halo, and high energy positrons are produced 
through leptonic and hadronic decays of the bosons. We investigate a possi- 
bility to detect the positron signal in upcoming experiments such as PAMELA 
and AMS-02. We found that the dark matter signal can be distinguished from 
the background in the PAMELA experiment when the dark matter mass is less 
than 120 GeV and the signal flux is enhanced due to a small scale clustering of 
dark matter. Furthermore, the signal from the dark matter annihilation can 
be detected in the AMS-02 experiment, even if such enhancement does not 
exist. We also discuss the invisible width of the Higgs boson in this model. 
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I Introduction 

The hierarchy problem in the Standard Model (SM) is expected to give a clue to 
explore the physics beyond the SM. This problem is essentially related to quadrat- 
ically divergent corrections to the Higgs boson mass, and we need a mechanism to 
avoid the divergences. To solve the problem, many scenarios have been proposed so 
far, for example supersymmetry, in which the divergences are completely removed. 
Other examples are scenarios with a low energy cutoff scale around a TeV such as 
Techni-color and TeV scale extra-dimension. 

The latter scenarios are, however, constrained by the electroweak precision mea- 
surements. From the analysis of higher dimensional operators at the cutoff scale, it 
has been found that the scale should be larger than roughly 5 TeV For such high 
energy cutoff, the hierarchy problem appears again: We still need the fine-tuning of 
a few percent level in the Higgs mass term in order to obtain the 100-200 GeV Higgs 
boson mass. This problem is called the little hierarchy problem. 

Recently the little Higgs model [21 [3] has been proposed for solving the little 
hierarchy problem. In this scenario, the Higgs boson is regarded as a pseudo Nambu- 
Goldstone boson. New particles such as heavy gauge bosons and a top-partner are 
introduced, and all quadratic divergences to the Higgs mass term completely vanish 
at one-loop level due to these particles' contributions. Thus, the fine-tuning of the 
Higgs boson mass is avoided even if the cutoff scale is around 10 TeV. 

The original little Higgs model is still strongly constrained by the electroweak 
precision measurements [1]. This is mainly due to the contributions to electroweak 
observables from new heavy gauge bosons, because their masses are much smaller 
than the cutoff scale. In particular, direct couplings among a new heavy gauge boson 
and SM particles give sizable contributions to the observables. As a result, masses 
of new particles have to be raised, and the fine-tuning of the Higgs boson mass is 
reintroduced. 

To resolve the problem, the implementation of the Z2 symmetry called T-parity 
to the model has been proposed [S]-[Z]. Under the parity, the new particles are 
assigned to be — charge (T-odd), while the SM particles have -|- charge (T-even). 
Thanks to the symmetry, dangerous interactions stated above are prohibited, and 
the masses of new particles can be lighter. 

Due to the T-parity, the lightest T-odd particle becomes stable and a good can- 
didate of the dark matter. This is an interesting feature of the model, because the 
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existence of the dark matter is now established by recent cosmological observations 
[S]. Since the hghtest T-odd particle is electrically and color neutral, and has a mass 
of O(IOO) GeV [5] in many little Higgs models with T-parity, these models provide a 
WIMP (weakly interacting massive particle) dark matter [9], and are able to account 
for the large scale structure of the present universe [TU] . 

In this paper, we study the dark matter phenomenology in the littlest Higgs model 
with T-parity [6|, [3, [11] . The relic abundance of the dark matter in the thermal relic 
scenario has already been evaluated [11], and it has been found that the mass of 
the dark matter consistent with the WMAP observation [S] is around a few hundred 
GeV. In this paper, we focus on the indirect detection of this dark matter using the 
cosmic positrons. The dark matter in the halo associated with our galaxy frequently 
annihilates and produces high energy particles, for example, positrons [I3]-[T6], anti- 
protons [T7], etc. Then high energy positron excess in the cosmic ray provides an 
opportunity to search for the dark matter signal. 

In the littlest Higgs model with T-parity, the dark matter candidate is a heavy 
photon, which annihilates mainly into weak gauge bosons [IT]. Positrons are pro- 
duced through decays of the bosons. Since the dark matter annihilation occurs in 
the s-wave and weak gauge bosons can produce high energy positrons through lep- 
tonic decays, the resultant positron flux is large and its spectrum becomes harder 
than that of background positrons originating in a secondary production from cosmic 
protons. This feature is quite different from the spectrum of a bino-like neutralino 
dark matter in supersymmetric models, which is expected to be much softer. In this 
paper we calculate the positron flux from the dark matter annihilation in the model, 
and estimate the possibihty to detect these positrons in future experiments such as 
PAMELA [IB] and AMS-02 US]. 

This paper is organized as follows. In the next section, we briefly review the 
littlest Higgs model with T-parity, in particular, focusing on the mass spectrum 
in the gauge-Higgs sector and interactions relevant to the calculation of the dark 
matter annihilation. We also present the thermal relic abundance of the dark matter. 
Calculation of the positron spectrum from the dark matter annihilation in the halo 
is performed in Sec. HI using a diffusion model. Results of the positron flux are 
shown in Sec. IV. We also present the x^-analysis in order to investigate a possibility 
to detect the positron signal in future experiments. Sec.V is devoted to summary 
and discussions including the Higgs decay into the dark matter. In Appendix, we 
consider constraints on the model from electroweak precision measurements, and 
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show that the entire region restricted by the WMAP observation can be consistent 
with the measurements. 



II Dark Matter in Littlest Higgs Model with T-parity 



We briefly review the littlest Higgs model with T-parity. In particular, we focus 
on the mass spectrum in the gauge-Higgs sector and interactions relevant to the dark 
matter in the model. We also calculate the relic abundance of the dark matter and 
present a parameter region of the model consistent with the WMAP observation [11] . 
For the general review of little Higgs models and some phenomenological aspects, 
see Refsf 



Littlest Higgs Model with T-parity 



The littlest Higgs model [3j is based on a non-linear sigma model describing an 
SU(5)/SO(5) symmetry breaking. The non-linear sigma fleld S is given as 



; 



where / is the vacuum expectation value associated with the symmetry breaking. 
The Nambu-Goldstone boson matrix, H, and the direction of the symmetry breaking 
in the non-linear fleld, Sq, are written as 
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An [SU(2) xU(l)]^ subgroup in the global symmetry SU(5) is gauged, which is broken 
down to the diagonal subgroup identifled with the SM gauge group (SU(2)LxU(l)y). 
Due to the presence of the gauge interactions (and Yukawa interactions if we intro- 
duce), the global symmetry SU(5) is not exact, and the particles in the H fleld 
become pseudo Nambu-Goldstone bosons. 

Fourteen (= 24 — 10) Nambu-Goldstone bosons are decomposed into represen- 
tations under the electroweak gauge group as Iq © 3o © "^±1/2 © 3±i. The flrst two 
representations are real, and become longitudinal components of gauge bosons when 
the [SU(2) xU(l)]^ is broken down to the SM gauge group. The representations 2±i/2 
and 3-1-1 are a complex doublet identifled with the SM Higgs fleld {H in Eq.([2])) and 
a complex triplet Higgs fleld ($ in Eq.([2])), respectively. 



3 



The kinetic term for the S field is given as 



Cy, = ^Tr 



where 



(3) 



(4) 



Here, Wj{Bj) are the SU(2)j(U(l)j) gauge fields and gj{g'j) are corresponding gauge 
coupling constants. The generators of the gauge symmetries Qj and Yj are 



la'' \ 
^ 0, 

\ y 

/ \ 

-^ , r2 = diag(2,2,2,-3,-3)/10, (5) 
^0 a"* y 

where cr" are the Pauli matrices. 

In terms of above fields, the symmetry of the T-parity [S]-IZ] is defined as the 
invariance of the Lagrangian under the transformation: 



yi = diag(3,3,-2,-2,-2)/10 
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where fl = diag(l, 1,-1, 1,1). As a result of the symmetry, the gauge coupling gi{g'i) 
must be equal to 5'2(fi'2), namely gi = 92 = V^g {g[ = g'2 = V^g'), where g{g') is 
nothing but the coupling constant of the SM SU(2)L(U(l)y) gauge symmetry. 

Since the Higgs boson is the pseudo Nambu-Goldstone boson, its potential is 
generated radiatively [31 [H] 



A 



V{H, $) = XfTi [$"f<l>] - fi^H^H + - (h^h) 



(7) 



Due to the little Higgs mechanism, quadratic divergent corrections do not contribute 
to the Higgs mass /x^ at 1-loop level, while the corrections do contribute to the triplet 
Higgs mass term. Main contributions to yU^ come from the logarithmic divergent 
corrections at 1-loop level and quadratic divergent corrections at 2-loop level. As 
a result, /i^ is expected to be smaller than while the triplet Higgs mass term 
is proportional to p . The quartet couphng A is determined by the 1-loop effective 
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potential from gauge and top sectors. Since both /i and A depend on parameters at 
the cutoff scale, we treat these as free parameters in this paper. 

We discuss the mass spectrum of gauge and Higgs bosons. This model contains 
four kinds of gauge fields W^, W2, Bi and B2 in the electroweak gauge sector. The 
combinations, W = (W^f + W2)/V^ and B = {Bi + i?2)/V2, correspond to the 
SM gauge bosons for the 811(2)^ and U(l)y symmetry. The other combinations, 
Wh = (Wi - W2) /V^ and Bh = {Bi - B2) / \pl, are additional gauge bosons, which 
acquire the masses of through the SU(5)/SO(5) symmetry breaking. After 

electroweak symmetry breaking, neutral components of Wf^ and B^ are mixed and 
form mass eigenstates An and Zn. The masses of the heavy bosons are obtained as 

1 



mz^ = -[^A + C + ^iA-Cy + AB^)^gf , 

mA„ = ^(^ + C-v'(^-C)2 + 4i?2) ^^/, (8) 

where A = g'^{P — f ^/4), B = gg'v'^/A and C = g'^{P/5 — f ^/4). The mixing angle 
between and Bh are given as 

25 gg'v^ 

tanOH = , — - - .rot 2 ^7^' 

A-C + J{A-Cf + AB^ 4:Pi9^ - g'V^) 



which is suppressed by 0{v/ f). In addition to these gauge fields, we have the triplet 
Higgs boson $ in this model, and its mass is given by m| = A/^ = 2m\p / v'^ , where 
irih is the mass of the SM Higgs boson and v (~ 246 GeV) is the vacuum expectation 
value of the Higgs field. New heavy gauge bosons and the triplet Higgs boson are 
T-odd particles, while SM particles are T-even. 

The mass spectrum of T-odd particles are determined by two parameters, the 
breaking scale (/) and the Higgs boson mass {mh). For instance, in the case of = 
120 GeV and / = 700 GeV, m^^ = 100 GeV, mw„iz„) = 450 GeV and m$ = 500 
GeV. As shown in Eq. ([8]) , the mass of the heavy photon is considerably lighter than 
other T-odd particles due to the small hypercharge. Thus its stability is guaranteed 
by the T-parity conservation and becomes a candidate of a non-baryonic cold dark 
matter. 

In addition to these new particles, top-partners are introduced in this model in 
order to cancel the quadratic divergent contribution to the Higgs mass term from the 
top quark loop diagrams. Due to the T-parity, three kinds of partners exist, namely 
T-even heavy top (T+) which is introduced for the cancellation, T-partners of heavy 
top and top quark (T_ and The mass spectrum of these particles depends not 



5 




yy{Z) Ah; 



W(Z) Ah' 



h Ah 



h A 




Figure 1: Feynman diagrams for the annihilation of the dark matter (Ah). 

only on / and rrih, but also on other model parameters. Since these particles do 
not play a significant role in the dark matter phenomenology, we do not discuss this 
sector here. 

Relic Abundance of Dark Matter 



The dark matter (Ah) in the model annihilates mainly into weak gauge bosons, 
W~^W~ , ZZ through the diagrams in which the Higgs boson propagates in the s- 
channel. The dark matter also annihilates into Higgs bosons if > m^. The 
Feynman diagrams for these processes are shown in Figl^. From Eqs.(l3]) and (171), 
interactions relevant to the annihilation are given as 



c { vh + 



4 2v 



(10) 



where c = —{gsmOu — g' cosOhY jA. We have used the unitary gauge H = (0,f + 
hY' /\/2. From the interactions, the annihilation cross sections of the dark matter 
turn out to be 
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^ There are also diagrams in which the T-partners of fermions are exchanged in the t-channel. 
These contributions are, however, negligible compared to those in FigUunless masses of T-partners 
are much smaller than 1 TeV. In fact, as discussed in Ref. [23, the cross sections for these processes 
are suppressed by masses of the T-partners, m^f^, and small hypercharges, Y'^ = (0.1)'*. 
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where v is the relative velocity between incident dark matters, and is the width 
of the SM Higgs boson. We take the non-relativistic limit (v — * 0) in the calculation, 
because the dark matter is almost at rest at the freeze-out temperature. 

The relic abundance of the dark matter is obtained by solving the following 
Boltzmann equation [23], 

where Y = n/s is the yield of the dark matter defined by the ratio of the dark matter 
density (n) to the entropy density of the universe 

and X = niAfj/T (T is the temperature of the universe). The Hubble parameter is 
given by H = l.QQgl^'^ m\^m-p\/ x"^ , where mpi = 1.22 x 10^^ GeV is the Planck mass. 
The yield in the equilibrium Fcq is written as 

Since the dark matter annihilates into the SM particles in the s-wave at the non- 
relativistic limit, the thermal averaged annihilation cross section ((crv)) is simply 
given by 

(crv) = (ty\^^ + aY\zz + (y^\hh > (14) 

After solving the Boltzmann equation, we obtain the present abundance of dark 
matter (Yqo). It is useful to express the relic density in terms of the ratio of the dark 
matter density to the critical density {VlK^ = mAgnh'^ / pc = mAHSoY^ooh"^ / pc), where 
= 1.1 X 10~5/i2 cm-3, h = O.ntom and Sq = 2900 cm'^. With a good accuracy, 
the solution of Eq. fll2l) is approximately given as 

= (15) 

^/^mpL(o-v) 

where Xf = niAjj/Tf is the freeze-out temperature for the dark matter and given 

1 /2 

as X/ = ln(X) — 0.51n(ln(X)) with X = 0.038 ■ {3/g* )'mpi^mAH{crv). Typically Xj 
takes a value, Xf ~ 23. 

The relic abundance of the dark matter in the thermal scenario is depicted in 
Figj2] as a contour map in the (/, m/i)-plane. The shaded thin area is the allowed 
region for the WMAP observation at 2a level, 0.094 < nh"^ < 0.129 [8]. The result 
obtained here is consistent with the previous calculation in Ref.[TT]. As shown in 
the figure, the breaking scale / is constrained to be less than about 2 TeV. 
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Figure 2: Contour plot of the thermal relic abundance for the dark matter {Qh'^) in the 
(/, m/i)-plane. The shaded thin area is the allowed region from the WMAP observation at 
2a level, 0.094 < Qh'^ < 0.129. 

The littlest Higgs model with T-parity is constrained from the electroweak preci- 
sion measurements. New physics contributions to the electroweak observables come 
from radiative corrections, because there is no tree-level effect due to the T-parity. 
As a result, the constraint to the model becomes weaker than that for the model 
without T-parity. The detailed analysis has been performed in Ref.[12]. We have 
repeated this analysis and present the result in Appendix. It is shown that the en- 
tire region is consistent with the electroweak precision measurements by choosing 
parameters in the top sector. 

Ill Propagation of Positron from Dark Matter Annihilation in Galaxy 

In the present universe, a dark matter makes up a halo associated with a galaxy. 
The distribution of the dark matter in the halo is given from a halo mass profile p(r) 
through the equation n{r) = p{f)/m. The profile is determined by observations of 
the rotational velocity of the galaxy and the motions of the dwarf galaxies with help 
of N-body simulations. Several models for the profile have been proposed [21] . For 
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our galaxy, we use the isothermal halo model in this paper, which is given as 

1 2/2 

= Po . t""/ 2 (GeV/cm=^) , (16) 
1 + r"^/r^ 

where r = |r| is the distance from the galactic center, po — 0.43 GeV/cm^ is the 
local halo density in the vicinity of the solar system, — 2.8 kpc is the core radius 
of the galaxy, and tq ~ 8.5 kpc is the distance between the galactic center and the 
solar system. 

In an indirect detection of dark matter, high energy particles from the dark 
matter annihilation are expected to be observed in the cosmic ray. Several kinds of 
annihilation products are produced, among which we focus on positrons [13]. We 
calculate the expected flux of the positrons at the earth from the annihilation. In 
evaluation of the flux, we need to take into account the propagation of positrons 
through the galaxy. We also address the background positrons originated from the 
secondary production of the cosmic ray. 

Production Rate of Positrons from Dark Matter Annihilation 

The dark matter annihilates mainly into weak gauge bosons and the Higgs bosons 
as discussed in the previous section. Positrons are produced through leptonic and 
hadronic cascade decays of these bosons. For W bosons, these processes are — > 
e'^u, fi^u ^ e^vvv or — > hadrons — > tt^ p^ e^. Decay branching 

ratio for the Higgs boson depends on its mass. When mh > 160 GeV, the Higgs 
boson mainly decays into weak gauge bosons, and positrons are produced by decays 
of weak bosons. 

Using the annihilation cross sections in Eg. (11 II) . the production rate of positrons 
from the annihilation is given as 

Q{E,f) = ln'ir) J: icrfv)(^] , (17) 

^ f=WW,ZZ,hh \ "'^ ) f 

where E is the energy of a positron and the coefficient 1/2 comes from the pair 
annihilation of the identical particles. The fragmentation function {dNf.+ / dE) f rep- 
resents the number of positrons with energy E produced from the final state /. The 
cascade processes for the positron production discussed above are encoded into these 
functions. 

The fragmentation functions are evaluated by a Monte-Carlo simulation such as 
the HERWIG code [25]. These functions for the weak gauge bosons can be parame- 
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terized by a single variable x = E/m. The fitting functions have been constructed in 
Ref . [2n] to reproduce results of the simulation. We use these functions in this paper. 
The fragmentation function for the Higgs boson is also obtained by assuming the 
dominance of gauge boson decay modes. 

Although we use the isothermal model in Eq. fll6l) . the high energy positron flux 
does not strongly depend on the choice of dark matter halo models. Main differ- 
ence among proposed models appears in the region around the galactic center, and 
positrons produced around this region can not reach the earth without the significant 
energy loss. 

Recently, the effect of inhomogeneity in the local dark matter distribution on the 
positron flux is discussed based on the N-body simulations. It is shown that the 
positron flux from the dark matter annihilation is enhanced if there are clumps of 
the dark matter in the vicinity of the solar system [27|. The effect is parameter- 
ized as a boost factor (BF), which is defined by the ratio of the signal fluxes with 
inhomogeneity and without inhomogeneity. 



where the region of the integration is taken to be ^ ~ (a few kpc)'^ around the solar 
system. The boost factor is larger than 1, and equal to 1 only if the density p is a 
constant. The value of the factor is expected to be in the range of 2 to 5 based on a 
hierarchical clustering scenario in the inflationary universe |28] . 

Positron Propagation in Galaxy 

Once positrons are produced in the dark matter annihilation, they travel in our 
galaxy under the influence of a tangled magnetic field. Since the typical strength 
of the magnetic field is a micro Gauss, the gyro-radius of the positron is much less 
than the galactic radius. Thus, the propagation can be treated as a random walk. 

We use a diffusion model for the propagation of positrons, in which the random 
walk is described by the following diffusion equation [151 , 



^JeAE,f) = K{E)V'U{E,r) + — [b{E)U{E,f)] + Q{E,r) , (19) 



where fe+{E,r) is the number density of positrons per unit energy, E is the energy 
of positron, K{E) is the diffusion constant, b{E) is the energy loss rate, and Q{E, r) 




(18) 
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is the source (positron injection) term discussed in the previous subsection. The flux 
of positrons with high energy {E ^ rrie) in the vicinity of the solar system is given 
from fe+{E, r) as 

^AE)=BF^U{E,r^) , (20) 

where Tq represents the coordinate of the solar system. 

There are two parameters in Eg. (1191) . One is the diffusion constant K{E) char- 
acterizing the tangled magnetic field of the galaxy. This parameter is evaluated by 
comparing the observations of the Boron to Carbon ratio in the cosmic ray with the 
result of simulations [29j. The parameter b{E) stands for the energy loss rate of 
positrons due to the inverse Compton scattering with cosmic microwave radiation 
(and infrared photons from stars) and synchrotron radiation with the magnetic field 
during the propagation in the galaxy [30]. This parameter is, therefore, determined 
by the photon density, the strength of the magnetic field and the Thomson scattering 
cross section. For both parameters, we use values adopted in Ref.[T5] 



K{E) = 3.3 X 10'' [3°-^ + (E/l GeVf^] (cm's^^) , 
b{E) = 10~^%E/1 GeVf (GeVs-^) . (21) 

The positrons from the dark matter annihilation are expected to be in the equi- 
librium in our galaxy, hence the number density fe+{E,r) is obtained by solving 
Eq. f|T9l) with the steady state condition dfe+/dt = 0. Furthermore, we impose the 
free escape boundary condition, namely the positron density drops to zero on the 
surface of the diffusion zone. This means that the positrons coming from the outside 
of the zone are negligible, while the positrons produced inside the zone contribute to 
the flux around the solar system due to the trapping by the tangled magnetic field 
|31j . The diffusion zone is usually assumed to be a cylinder with the half- height (L) 
and radius (R), which are set to be L = 4 kpc and R = 20 kpc in this paper. 

The high energy positron flux at the earth does not strongly depend on the choice 
of the diffusion zone, because the positrons we observe are produced within a few 
kpc around the solar system. In fact, the distance in which positrons travel without 



significant energy loss is estimated as r ~ ^EK{E)/b{E) ~ 0.74 x (E/lOO GeV)~°-' 
kpc. Positrons far from the earth lose their energies during the propagation, and 
consequently they only contribute to the low-energy part of the flux. 

The flux obtained from Eq. (!20|) does not correspond exactly to the one observed 
on the top of atmosphere. The flux is modified due to interaction with the solar 
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wind and the magneto-sphere [32]. However, the modulation effect is not important 
when the energy of a positron is above 10 GeV. Furthermore, the effect is highly 
suppressed in the positron fraction, which is defined by a ratio of the positron flux 
to the sum of positron and electron fluxes, i.e. ^e+/{^e+ + ^e-)- 

Background Positrons 

In the detection of the positrons from the dark matter annihilation, the main 
background comes from high energy positrons in the cosmic ray. High energy 
positrons are produced as secondary particles in the collision between hydrogen and 
helium in interstellar medium and primary particles in cosmic ray accelerated by 
the shock wave in supernovas. The flux of these positrons are obtained by simula- 
tions, in which a diffusion model is also used. The result of simulations agrees with 
measurements of the low-energy positron flux in the cosmic ray [33]. The fitting 
functions for high energy positrons, primary electrons, and secondary electrons have 
been constructed [15] , 

- .^i^^L^.. (OeV-cm-s-sr-) , 

- l + 110i.^.^r600r^.^ + 580i.^.^ (GeV-cm-s-sr-) , 

^'fr\E) = (GeV-^cm-^s-isr-i) , (22) 

^+ ^ ' 1 + 650^2-3 + 1500^4-2 ^ ^ ' ^ ' 

where E is in unit of GeV. The first one, ^^P""'^ is the flux of the primary electrons, 
while the second and third ones, ^'■"^"^ and are the secondary electron and 

positron fluxes, respectively. 

IV Positron Signal from Dark Matter Annihilation in Halo 

We are now in position to discuss the positron signal from the dark matter 
annihilation. The signal positron flux is evaluated in Eg. (1201) . while expected positron 
and electron background are given in Eqs. (!22|) . In order to show how the dark matter 
annihilation can modify the positron energy spectrum in the cosmic ray, we have 
chosen seven sample-points (I to VII) in the parameter space of the little Higgs model 
with T-parity as in FigjSl Parameters / and rrih and masses of various particles are 
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Figure 3: Sample-points for depicting the positron fraction from the dark matter anni- 
hilation. The shaded thin area is the allowed region from the WMAP observation at 2a 
level. In the table below the figure, the details about model parameters in each point are 
shown. 

listed in the table below the figure. All points satisfy the WMAP condition, namely 
the present dark matter abundance is explained by the thermal relic scenario. As 
seen in the figure, there are two branches: the upper branch (U-branch) and lower 
branch (L-branch). In the U-branch, the Higgs boson mass is larger than twice the 
dark matter mass, rrih > 2171 a^, while rrih < 2m in the L-branch. 

In FigHl the positron fraction, $e+ / ('^'e+ + '^'e-); is shown as a function of positron 
energy. In the left figure, the results for the points in the U-branch (I to III) are 
depicted, while those in the L-branch (IV to VI) are in the right figure. The point VII 
can be regarded as a sample on both U- and L-branches, and its result are shown in 
both figures. We used the boost factor BF = 5. The expected background positron 
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Figure 4: The positron fraction as a function of positron energy E. For comparison, the 
expected background fraction is also shown in these figures. In the left figure, the fraction 
in the U-branch (I to III and VII) are depicted, while those in the L-branch (IV to VI and 
VII) are in the right figure. In both figures, the boost factor BF = 5 is used. 



fraction is also shown for comparison. 

The upcoming experiments such as PAMELA and AMS-02 have good sensitivities 
in a broad range for a positron energy 10 GeV < E < 270 GeV. In order to discuss 
the possibility for detection of the dark matter signal in the future experiments, we 
perform the x^-analysis developed in Ref.[T6]. For this purpose, we need to know 
the expected signal and background events in future experiments for each parameter 
point of the model (/, rrih). The is defined as 

X^ = E ^ ' ,W ^ > (23) 

i jVj 

where the sum is taken over energy bins, Nf'^^'' is the number of positron events 
observed in the i-th bin and Nl is the number of events expected from the back- 



ground contribution in the bin. Following the Ref . [TB] , we chose 22 bins in the range 
between 10 GeV < E < 270 GeV, 

A[logio(^/lGeV)] = 0.06 , {E < 40GeV) , 
A[logio(^/lGeV)] = 0.066 , {E > 40GeV) . (24) 

In our analysis, we use the acceptance of PAMELA and AMS-02 to be 20.5cm^sr 
and 450cm^sr, respectively, assuming three years of data-taking. 

In Figini the contour plot of the is depicted in (/, m/i)-plane. The left figure 
is the ill the PAMELA experiment with BF = 5, while the right one is the 



14 



mAH(GeV) mAH(GeV) 




f (TeV) f (TeV) 

Figure 5: The contour plot of the in (/, ?TT-/j)-plane. The left figure is the ^ in the 
PAMELA experiment with BF = 5, while the right one is the in the AMS-02 with 
BF = 2. For comparison, the constraint from the WMAP observation is also shown 
as a shaded region. The values of x^, 30.8, 33.9 and 40.3, correspond to the statistical 
significance for the detection of the signal at the 90%, 95% and 99% confidence levels. 



in the AMS-02 with BF = 2q The statistical significance for the detection of the 
signal at the 90%, 95% and 99% confidence levels correspond to the values of x^, 
30.8, 33.9 and 40.3. In these figures, the constraint from the WMAP observation is 
also shown as a shaded region. Inside the region, predicted thermal relics of the dark 
matter are smaller, while the relics outside the region are larger than that observed 
in the WMAP observation. If we like to consider parameter region inconsistent with 
the WMAP observation, deviation from the thermal relic scenario is necessary such 
as a non-thermal production of dark matter or entropy production at late time. 

The strange behavior around / > 1 TeV and rrih ~ 200-400 GeV in the figures 
is due to the annihilation mode into two Higgs bosons. When this annihilation 
channel is opened {rriAu > Trih) and the Higgs boson mass is larger than twice the W 
boson mass {rriH > 2m\y), high energy positrons are produced through the process 
AhAh —>■ hh ^ WWWW . These positrons have a hard spectrum, and enhance 
the possibility to detect the dark matter signal. On the other hand, when the Higgs 
boson mass is less than that of two W bosons {mfj < 2miy), the effect becomes 
negligible. Since positrons are produced through cascade decays of b-quarks, the 
resultant spectrum in the process becomes very soft. 

^The value of is proportional to BF^, thus the extension of the result in other values of the 
boost factor is straightforward. 
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Figure 6: plot along the center value {Qh'^ = 0.112) of U- (left figure) and L-branch 
(right figure) as a function of the Higgs boson mass nih and the dark matter mass ttiah (see 
the top axis). In both figures, the of PAMELA with BF = 5, AMS-02 with BF = 2, 
and BF = 1 are depicted. 



From the left figure, we see that the dark matter signal is clearly distinguished 
from the background in the PAMELA experiment when the breaking scale / is less 
than 1 TeV for BF = 5. On the other hand, from the right figure, almost all 
interesting area including the WMAP region is covered in AMS-02 with BF = 2. 
Furthermore, the value of in the AMS-02 with BF = 1 is quite similar to the plot 
in the PAMELA with BF = 5. Thus it is possible to detect the signal in AMS-02 
even if there is no enhancement from the boost factor. 

The plot along the center value {Qh"^ = 0.112) of U- and L-branch as a function 
of the Higgs boson mass m/j and dark matter mass rriA^ are presented in FiglHl The 
result of the U-branch case is shown in the left figure, while that of the L-branch is 
in the right figure. In both figures, the of the PAMELA with BF = 5, AMS-02 
with BF = 2 and BF = 1 are depicted. For the reference, the line = 33.9 
(corresponds to the 95% confidence level) is shown. The decreasing behavior of 
along with increasing rrih is due to the fact that a number density of the dark matter 
is decreasing as niAfj is increasing. If the boost factor is around 5, the PAMELA 
experiment has a potential to detect the dark matter signal when nih < 300 GeV 
{mAfj < 120 GeV) in the U-branch case or rrih < 150 GeV {rriAH < 120 GeV) in the 
L-branch. Furthermore, these regions can be covered in AMS-02, even if BF = 1. 

Finally, we show the 95% confidence level contour within the WMAP constraint 
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Figure 7: 95% confidence level for the statistical significance of the dark matter 
detection within the WMAP constraint in PAMELA and AMS-02 experiments. The 
region above the line can be distinguished from the background spectrum. The plot 
is depicted as a function of the breaking scale / (or the dark matter mass m^^ on 
the top axis) and the boost factor BF. The lines for U- and L-branches coincide. 



in (/ or myi^,i?F)-plane. The region above the line can be distinguished from the 
background in each experiment. Although both results in U- and L-branch cases 
are depicted in the figure, the contour lines are almost degenerate in this parameter 
space. In both branches, the annihilation modes are completely dominated by WW 
and ZZ bosons, and the positron production cross sections are the same once the 
WMAP constraint is applied. From the figure, we see that the signal may be detected 
in the PAMELA experiment when / < 830 GeV (m^^ < 120 GeV) and BF > 5. 
On the other hand, the AMS-02 experiment will cover a wide range of the parameter 
space including the region with BF = 1. 

Here, we comment on the positron excess recently reported by HEAT collabora- 
tion [34J. In the measurement, the excess of high energy positrons (1 GeV < E < 
30 GeV) has been observed. If the excess is due to the dark matter annihilation, 
its annihilation cross section should be large (crv ~ lO^^^cm^sec"^) unless the boost 
factor is large {BF ~ 50-100). Such large annihilation cross section is difficult to 
satisfy the WMAP constrainlo. It is therefore unlikely that the excess observed at 



^The cross section (crv ^ 2 x 10 ^^cm'^sec ^) is required in order to explain WMAP result in 
the thermal relic scenario 
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Figure 8: Branching ratio of the Higgs decay process into two dark matters along 
the center value of the U-branch as a function of the Higgs boson mass. 



the HEAT experiment is explained by the signal of the dark matter in this mode^ 



V Summary and Discussions 

We have studied the possibility to detect the dark matter in the littlest Higgs 
model with T-parity in future cosmic positron experiments. High energy positrons 
are produced from the dark matter annihilation through weak gauge boson decays. 
The resultant positron spectrum becomes hard, and the indirect detection of the dark 
matter is promising. We have performed the x^-analysis to evaluate a confidence level 
to detect the dark matter signal in upcoming experiments. We have found that the 
signal will be detected in the PAMELA experiment, when the dark matter mass is 
less than 120 GeV and the boost factor is around 5 within the WMAP constraint. 
The region itiah < 120 GeV (/ < 830 GeV) corresponds to nih < 300 GeV in the 
U-branch, and rrih < 150 GeV in the L-branch. In the AMS-02 measurement, these 
regions can be covered even if there is no enhancement from the boost factor. 

The positron spectrum in this model has a different feature compared to a bino- 
like neutralino dark matter in supersymmetric models. In the supersymmetric case, 
positrons are mostly produced from bottom quark decays, so that its spectrum be- 
comes softer than the dark matter in the present model. Furthermore, the positron 

^There are large uncertainties on the positron flux due to a solar modulation at GeV energy, 
which may be responsible to the HEAT excess [32] . 
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production cross section is much smaller, because the bino-like dark matter anni- 
hilates in the p-wave. Therefore, the indirect detection of the dark matter in the 
littlest Higgs model with T-parity is easier than the supersymmetric case. 

Finally, we discuss possible effects of the dark matter on the Higgs phenomenol- 
ogy. In the U-branch, the decay of a Higgs boson into a pair of the heavy photons 
contributes to the invisible width. The branching fraction of the process is shown in 
Figj8]as a function of the Higgs boson mass. We see that the ratio is at a few percent 
level. Although this value seems to be beyond the currently estimated sensitivity 
to the total width measurement of the Higgs boson in LHC [33] and ILC [SB], the 
measurement might be possible at a future muon collider [27]. 

We have concentrated on the positron signal of the dark matter in this paper. 
It is also interesting to consider other ways to search for the dark matter signals 
such as direct detection, and indirect detections using neutrinos, gamma-rays and 
anti-protons. 
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Appendix Constraints from Electroweak Precision Measurements 

In Appendix, we consider constraints on the Littlest Higgs model with T-parity 
from electroweak precision measurements. We follow the procedure in Ref. [T2] using 
S, T and U parameters [38j. In that paper, it has been shown that main contributions 
to these parameters come from the top-sector in the model and custodial-symmetry 
violating effects from heavy gauge boson loops. 

Contribution from the top-sector depends not only on the parameter / but also 
on a new parameter, R = A1/A2. Using Ai and A2, masses of the top quark and 
its T-even partner are given as rrit = AiA2f/(A^ + X^Y^"^ and = (A^ + X^Y^'^f. 
Instead of Ai, A2, and /, we take m^, R, and / as free parameters. We have calculated 
contributions from the top sector to the S, T and U parameters, and confirm the 
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results in Ref. [12], which are given by 
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Tr 



Ur = ^ 



-7T +log 
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mt 




(25) 



We have also calculated contributions from heavy gauge boson loops in the gen- 
eral i?^ gauge. We have obtained ^ independent results. Final expression, however, 
is not consistent with the result in Ref. [T2j, in which the Landau gauge has been 
used. Detail of our calculation is presented in the foUowings. 

In order to define the gauge, we first consider the mixing term between gauge 
bosons and derivatives of Nambu-Goldstone (NG) bosons up to appropriate order of 
v/f. The non-linear sigma field S is expanded by the NG fields as 

S = exp[22((n)+5n)//]So 

= (constant) + (terms oc 611) + (terms oc 6lP) + ■ ■ ■ , 
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We used the notation adopted in Ref. [12] for the NG fields in the 611 matrix. It is 
well known that the constant term is obtained in an exact manner as 



(constant) 
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where s = sm(\^v/ f), c = cos{y/2v/ f). Thanks to the Feynman formula [39], we 
can also derive terms proportional to 611 exactly as follows, 



(terms oc 611) 



(29) 



The integration by the parameter a can be exactly performed. Due to the complex 
form of the result, we omit to write it down. Using these exact expressions in Eqs.( 
and (l29l) . the mixing term induced from the kinetic term in Eq.Q is written as 



2v 
8 



f V2v/ 
' 1 sc 
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,(30) 



where Wh = (Wjj + iW^)/V2, and Zh = Wfj. We neglect the effect of U(l) gauge 
interactions for simplicity {g' = 0). 

Due to the electroweak symmetry breaking, the NG mode absorbed in the longi- 
tudinal component of the heavy gauge boson Wh (Zh) is given by the combination 



of u and 



AO 



p, and rj). Thus, "would-be NG" modes are written as 
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Note that these modes become vr", uj, and uj^ at the leading order of f //. Normaliza- 
tion constants, A^^. {x = vr, uj, and uj^), are determined by considering kinetic terms 
of NG bosons. The kinetic terms are also obtained exactly as 
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As seen from the equation, we find that kinetic terms of tt, uj, (p, uj^, 0p, rj, and 0° are 
not canonically normalized due to the electroweak symmetry breaking. Therefore, 
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we have to redefine these NG fields by normahzed ones. After the redefinition, we 
found the normalization constants, A^^;, in Eg. (1311) are given as 
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(33) 



Once we obtain canonically normalized "would-be NG" modes, we can determine 



gauge fixing functions to cancel the mixing term in Eg. (1301) . Those functions are 

(34) 
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Using these functions, mass terms of the "would-be NG" modes turn out to be 



= ^rnl {n-f + ^ml^Co^Co + | (m^^ + AM') {u')' , (35) 

NG mass 

where = g'p{l - c)/4, m^r^ = g'p{3 + c)/4, and AM' = g'f\l - c)V8. 
Hence, with ^ = 1, "would-be NG" masses coincide with those of corresponding 
gauge bosons at any order of vj f. 

We are now in position to calculate the contribution to the T parameter from 
heavy gauge boson loops. With expressions of "would-be NG" modes in Eg. (l?Tl) . 
interactions relevant to the calculation are given as follows. 
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In order to calculate the coefficients, C, in Eg. (1401) . we need not only the linear terms 
in Eg. (12^ . but also those proportional to 5112. gj^(.g jg difficult to obtain the 5112 
terms in an exact manner, we have calculated these terms with the expansion oivjj. 
Up to the order of {vj f)^, the coefficients turn out to be 
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(42) 



Interactions in >CvlVhVh '^VlVlVhVh gauge self-interactions coming from 
gauge kinetic terms, while ghost interactions in /^chost are derived from the gauge 
fixing functions in Eg. fl51|) . Interactions in >CvlNgnG) -^VlVlNgng and ^VlVhNG have 
been obtained by expanding the kinetic term in Eg. ([3]) in terms of NG fields. 
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Figure 9: Heavy gauge boson loop diagrams contributing to the T parameter. 



There are six kinds of loop diagrams involving heavy gauge and NG bosons con- 
tributing to the T parameter as shown in FigJH Those are diagrams with three points 
gauge self-interactions (VlVrVh), four points gauge self-interactions (VlVlVrVh), 
ghost interactions (Ghost), three points NG gauge interactions (VlNGNG), four 
points NG gauge interactions (Vl VlNGNG), and interactions in ^VlVhNG (VlVrNG). 
In all diagrams except the last one, the contributions to the T parameter come from 
the mass difference between Wh and Zh- In the last diagram, difference among 
coefficients in Eq.( l4Tl) is the source of the T parameter. 

The leading logarithmic divergent contribution from each diagram is 
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After summing all contributions, we have found that the logarithmic divergent cor- 
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Figure 10: Contributions to S and T parameters in the Littlest Higgs model with T-parity 
(solid lines) for < i? < 2 and / along the center value of the U-branch (left figure) and 
L-branch (right figure). From bottom to top, lines correspond to i? = 0, 0.5, 1, 1.5 and 
2. From left to right, / corresponds to the point I, II, III, and VII for U-branch, IV, V, 
VI, and VII for L-branch in FigjSl Constraints on S and T parameters from electroweak 
precision measurements are also shown at 68% and 99% confidence level. The plot assumes 
U = 0. 



rection is completely canceled, 

This result differs from Eq.(3.6) in Ref.[T2]. Therefore, the contribution to the T 
parameter at this order of v/ f comes from a finite term. 

In addition to the loop contribution, we expect that the following operator can 
arise at the cutoff scale [T2l. 



which gives the contribution to the T parameter at the same order of the finite term, 

nCutoff) ^ (^) . (46) 

The value of Sc is determined by the UV-completion of the Littlest Higgs model with 
T-parity, and it is naturally expected to be 0{1). By comparing the contribution 
from the top sector with those from heavy gauge boson loops and the cutoff scale 
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Figure 11: Constraints on the model parameters, R and /, at 68% and 99% confidence 
level. At each point in these figures, the Higgs mass is determined to satisfy the WMAP 
constraint on the U-branch (left figure) and L-branch (right figure). 



operator, the latter contributions turns out to be small as long as the finite term in 
Eq. fH^ does not have a large coefficient. In fact, T( Cutoff) in Eq. fH6l) is 8% of Tt 
in Eg. (1251) . when we take / = 1 TeV, i? = 1, and 5c = 1. Therefore, we neglect these 
contributions to derive electroweak precision constraints on the model. 

The contributions to S and T parameters from the Littlest Higgs model with T- 
parity are depicted in FiglTOl Constraints on S and T parameters from electroweak 
precision measurements at 68% and 99% confidence level are also shown. Solid 
lines are predictions of the model for < i? < 2 and / along the center value 
of the U-branch (left figure) and L-branch (right figure). For depicting contours 
of the constraints from the measurements, we have used three experimental values 
following papers [ID], W boson mass {rrtw = 80.412 ± 0.042 GeV), weak mixing 
angle [sm^ 9^^fj = 0.23153 ± 0.00016), and leptonic width of the Z boson (F; = 
83.985 ±0.086 MeV) [41j. We have also used values of the fine structure constant on 
the Z pole {a-^{mz) = 128.950 ± 0.048) and the top quark mass (m^ = 172.7 ± 2.9 
GeV) [42J. The origin of the S-T plane is fixed by using the reference Higgs mass, 
rrih = 100 GeV. In this figure, we can see that, even for a large value of /, there 
is a parameter region consistent with the electroweak precision measurements with 
appropriate choice of R. In the case of large /, the Higgs boson becomes heavy and 
gives large contributions to the S and T parameters. These contributions, however, 
can be canceled by those from the top loops {t and T+) due to the large mixing angle 
iR>l). 

In FigJTT| constraints for the model parameters, R and /, are shown at 68% and 
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99% confidence level. At each point in these figures, the Higgs mass is determined 
to satisfy the WMAP constraint on the U-branch (left figure) and L-branch (right 
figure). Possible range of R can be determined by requiring Ai and A2 couplings 
to be within the perturbative range, which results in 0.2 ^ R^5. As seen in these 
figures, the entire allowed region can be consistent with the electroweak precision 
measurements. 



References 

[1] R. Barbieri and A. Strumia, Phys. Lett. B 433 (1998) 63; R. Barbieri and 
A. Strumia, |arXiv:hep-ph70007265] 



[2] N. Arkani-Hamed, A. G. Cohen and H. Georgi, Phys. Lett. B 513 (2001) 
232; N. Arkani-Hamed, A. G. Cohen, E. Katz, A. E. Nelson, T. Gregoire and 
J. G. Wacker, JHEP 0208 (2002) 021. 

[3] N. Arkani-Hamed, A. G. Cohen, E. Katz and A. E. Nelson, JHEP 0207 (2002) 
034. 

[4] C. Csaki, J. Hubisz, G. D. Kribs, P. Meade and J. Terning, Phys. Rev. D 67 
(2003) 115002; J. L. Hewett, F. J. Petriello and T. G. Rizzo, JHEP 0310 (2003) 
062; C. Csaki, J. Hubisz, G. D. Kribs, P. Meade and J. Terning, Phys. Rev. D 

68 (2003) 035009; T. Gregoire, D. R. Smith and J. G. Wacker, Phys. Rev. D 

69 (2004) 115008; M. C. Chen and S. Dawson, Phys. Rev. D 70 (2004) 015003; 
Z. Han and W. Skiba, Phys. Rev. D 72 (2005) 035005; W. Kilian and J. Renter, 
Phys. Rev. D 70 (2004) 015004. 

[5] H. C. Cheng and I. Low, JHEP 0309 (2003) 051. 

[6] H. C. Cheng and I. Low, JHEP 0408 (2004) 061. 

[7] I. Low, JHEP 0410 (2004) 067. 

[8] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148 (2003) 
175; C. L. Bennett et al, Astrophys. J. Suppl. 148 (2003) 1. 

[9] For reviews, 

G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267 (1996) 195; 
L. Bergstrom, Rept. Prog. Phys. 63, (2000) 793; G. Bertone, D. Hooper and 
J. Silk, Phys. Rept. 405 (2005) 279; C. Munoz, Int. J. Mod. Phys. A 19 (2004) 
3093. 



27 



[10] J. R. Primack, Nucl. Pliys. Proc. Suppl. 124 (2003) 3. 

[11] J. Hubisz and P. Meade, Phys. Rev. D 71 (2005) 035016, (For the correct 
paramter region consistent with the WMAP observation, see the figure in the 
revised vergion, |hep-ph/0411264 V3). 

[12] J. Hubisz, P. Meade, A. Noble and M. Perelstein, JHEP 0601 (2006) 135. 

[13] S. Rudaz and F. W. Stecker, Astrophys. J. 325 (1988) 16; J. R. Elhs, R. A. Flo- 
res, K. Freese, S. Ritz, D. Seckel and J. Silk, Phys. Lett. B 214 (1988) 403; 
A. J. Tylka, Phys. Rev. Lett., 63 (1989) 840; M. S. Turner and F. Wilczek, 
Phys. Rev. D, 42 (1990) 1001; A. J. Tylka, Phys. Rev. D 43 (1991) 1774. 



[14] E. A. Baltz, J. Edsjo, K. Freese and P. Gondolo, arXiv: astro-p h/02 11239 
G. L. Kane, L. T. Wang and T. T. Wang, Phys. Lett. B 536 (2002) 263; W. de 
Boer, C. Sander, M. Horn and D. Kazakov, Nucl. Phys. Proc. Suppl. 113 (2002) 
221; G. L. Kane, L. T. Wang and J. D. Wells, Phys. Rev. D 65 (2002) 057701; 

E. A. Bahz, J. Edsjo, K. Freese and P. Gondolo, Phys. Rev. D 65 (2002) 063511; 
D. Hooper and G. D. Kribs, Phys. Rev. D 70 (2004) 115004; S. Profumo and 
P. Ullio, JCAP 0407 (2004) 006. 

[15] E. A. Baltz and J. Edsjo, Phys. Rev. D 59 (1999) 023511. 

[16] D. Hooper and J. Silk, Phys. Rev. D 71 (2005) 083503. 

[17] J. Silk and M. Srednicki, Phys. Rev. Lett. 53 (1984) 624; S. Rudaz and 

F. W. Stecker, Astrophys. J. 325 (1988) 16; A. Bottino, F. Donato, N. Fornengo 
and P. Salati, Phys. Rev. D 58 (1998) 123503; L. Bergstrom, J. Edsjo and P. Ul- 



lio, arXiv:astro-ph/9906034; F. Donato, N. Fornengo, D. Maurin, P. Salati and 
R. Taillet, Phys. Rev. D 69 (2004) 063501. 

[18] M. Circella [PAMELA Collaboration], Nucl. lustrum. Meth. A 518 (2004) 153; 
P. Spillantini, Nucl. lustrum. Meth. B 214 (2004) 116. 

[19] F. Barao [AMS-02 Collaboration], Nucl. lustrum. Meth. A 535 (2004) 134. 

[20] M. Schmaltz and D. Tucker-Smith, larXiv:hep-ph/0502182^ M. Perelstein, 



arXiv:hep-ph/0512128 



[21] G. Burdman, M. Perelstein and A. Pierce, Phys. Rev. Lett. 90 (2003) 241802 
[Erratum-ibid. 92 (2004) 049903]; T. Han, H. E. Logan, B. McElrath and 
L. T. Wang, Phys. Rev. D 67 (2003) 095004; M. Perelstein, M. E. Peskin and 
A. Pierce, Phys. Rev. D 69 (2004) 075002. 



28 



[22] A. Birkedal, A. Noble, M. Perelstein and A. Spray, |arXiv:hep-ph/0603077 



[23] E. W. Kolb and M. S. Turner, The Early Universe, (Addison- Wesley, Reading, 
MA, 1990). 

[24] J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 462 (1996) 563; 
J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 490 (1997) 493; 
B. Moore, S. Ghigna, F. Governato, G. Lake, T. Quinn, J. Stadel and P. Tozzi, 
Astrophys. J. 524 (1999) L19; J. Diemand, B. Moore and J. Stadel, Mon. Not. 
Roy. Astron. Soc. 353 (2004) 624. 

[25] G. Corcella et ai, interfering gluons (including supersymmetric processes)," 
JHEP 0101 (2001) 010. 

[26] J. Hisano, S. Matsumoto, O. Saito and M. Senami, |arXiv:hep-ph/0511118l 

[27] J. Silk and A. Stebbins, Astrophys. J. 411 (1993) 439; L. Bergstrom, J. Edsjo 
and P. Gondolo, Phys. Rev. D 59 (1999) 043506. 

[28] V. Berezinsky, V. Dokuchaev and Y. Eroshenko, Phys. Rev. D 68 (2003) 103003. 

[29] A. W. Strong and I. V. Moskalenko, Astrophys. J. 509 (1998) 212; 
I. V. Moskalenko and A. W. Strong, Phys. Rev. D 60 (1999) 063003; D. Maurin, 
F. Donato, R. Taillet and P. Salati, Astrophys. J. 555 (2001) 585; D. Maurin, 
R. Taillet and F. Donato, Astron. Astrophys. 394 (2002) 1039; I. V. Moskalenko, 
A. W. Strong, S. G. Mashnik and J. F. Ormes, Astrophys. J. 586 (2003) 1050; 
D. Hooper, J. E. Taylor and J. Silk, Phys. Rev. D 69 (2004) 103509. 

[30] M. S. Longair, High-energy astrophysics. Vol. 2 Cambridge University Press, 
New York, 1994. 

[31] A. Barrau, G. Boudoul, F. Donato, D. Maurin, P. Salati and R. Taillet, Astron. 
Astrophys. 388 (2002) 676. 

[32] L. J. Gleeson and W. I. Axford, Astrophys. J. 149 (1967) L115. 

[33] I. V. Moskalenko and A. W. Strong, Astrophys. J. 493 (1998) 694. 

[34] S. W. Barwick et al. [HEAT Collaboration], Phys. Rev. Lett. 75 (1995) 390. 
S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482 (1997) L191; 

[35] V. DroUinger and A. Sopczak, Eur. Phys. J. directC 3 (2001) Nl. 

[36] N. Meyer and K. Desch, Eur. Phys. J. C 35 (2004) 171. 

[37] M. M. Alsharoa et al. [Muon Collider/Neutrino Factory Collaboration], Phys. 
Rev. ST Accel. Beams 6 (2003) 081001. 



29 



[38] M. E. Peskin and T. Takeuchi, Phys. Rev. D 46 (1992) 381. 

[39] See, for exapmle, p. 12, Yuri Makeenko, Method of Contemporary Gauge The- 
ory" (Cambridge university press). 



[40] W. J. Marciano, |arXiv:hep-ph/0003181t M. Perelstein, M. E. Peskin and 
A. Pierce, Phys. Rev. D 69 (2004) 075002. 

[41] [ALEPH Collaboration], Phys. Rept. 427 (2006) 257. 

[42] [CDF Collaboration], [i^iv:hep-ex/050709l| 



30 



